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1.0  Introduction 


Multiple-Input,  Multiple-Output  (MIMO)  radar  concepts  have  recently  been  proposed  for  both 
airborne  and  surface-based  radar  applications  [1-13].  In  a  MIMO  radar,  Nt  transmit  antennas  each  emit 
a  unique  waveform.  These  waveforms  propagate  through  the  environment,  where  they  are  reflected  by 
targets  and  other  objects,  and  subsequently  received  by  Nr  receive  antennas.  Each  receiver  is  connected 
to  a  bank  of  filters,  with  one  filter  tuned  to  each  of  the  transmitted  waveforms.  Collectively,  these  filter 
outputs  form  N,  ■  Nr  data  channels  that  can  then  be  combined  as  desired,  e.g.,  to  maximize  Signal  to 
Interference  plus  Noise  Ratio  (SINR). 

In  the  literature,  MIMO  radars  have  been  shown  to  provide  a  number  of  potential  benefits  as 
compared  with  normal  phased-array  radars  (i.e.,  radars  that  emit  a  single  coherent  transmit  beam  that  is 
later  received  using  [l]an  array  of  receive  antennas).  These  benefits  include:  increased  angle  resolution 
[1-6],  increased  Doppler  resolution  [1,7],  reduced  clutter  levels  (and  related  hardware  requirements)  [1], 
sharper  airborne  radar  clutter  notches  [2,  7,  1 3],  and  Lower  Probability  of  Intercept  [1]. 

To  achieve  these  benefits,  most  of  the  literature  has  assumed  the  use  of  “orthogonal”  waveforms, 
hi  practice,  however,  such  waveforms  do  not  exist  and  must  be  approximated.  This,  in  turn,  has 
prompted  research  into  the  area  of  waveform  design  [14-16],  with  such  research  producing  a  variety  of 
frequency  offset  waveform  concepts.  Among  these,  Doppler  Division  Multiple  Access  (DDMA),  in 
which  waveforms  are  offset  in  frequency  within  the  Doppler  band  of  the  radar,  offers  particular  promise. 

After  describing  the  MIMO  radar  signal  environment  (Section  2.0),  we  briefly  review  the  theory 
underlying  Doppler-offset  waveforms  (Section  3.0).  In  particular,  we  describe  how  DDMA  waveforms 
achieve  “quasi-orthogonality”  by  exploiting  prior  knowledge  of  the  signal  environment.  Then,  in  Section 
3.1,  we  discuss  the  range  and  Doppler  ambiguities  associated  with  DDMA  waveforms,  hi  Section  3.2, 
these  ambiguities  are  shown  to  lead  to  “phantom  blind  speeds”  and  other  problems.  Section  3.3  analyzes 
these  ambiguities,  and  shows  how  they  are  related  to  array  topology  and  waveform  parameters.  Then,  in 
Section  4.0,  we  describe  a  modified  DDMA  waveform  approach  that  enhances  both  target  detection  and 
clutter  cancellation  in  airborne  MTI  systems  and  related  radar  applications. 

2.0  MIMO  Radar  Principles 

Consider  an  array  of  IV,  transmit  antennas  and  Nr  receive  antennas,  hi  MIMO  radar,  each 
transmit  antenna  radiates  an  independent  waveform  s„ .  Here,  s„  denotes  an  Lx  1  vector  containing 
complex  baseband  samples  of  the  «th  transmit  waveform.  Collectively,  these  /V,  transmit  waveforms 
form  an  N,xL  signal  matrix,  S  =  [s,  s2  •••  sN  ]  .  These  N:  waveforms  propagate  through  the 
environment,  where  they  are  reflected  by  both  targets  and  clutter.  In  particular,  the  waveforms  reflected 
from  point-target  t  will  create  a  signal  that  is  proportional  to  a^(0,  )S,  where  (0, )  denotes  the 
N ,  x  1  response  vector  for  the  transmit  array.  Here,  0,  represents  a  set  of  target  signal  parameters 
describing  the  target  t.  In  typical  airborne  radar  systems,  for  example,  0,  =  {$,,£}  where  6t  and  £  are 
the  azimuth  and  depression  angles  from  the  radar  to  target  t,  respectively. 

The  reflected  signal  then  propagates  through  the  environment  and  is  received  by  an  array  of 
antennas.  Assuming  each  receiver  collects  R»L  samples,  the  NrxR  sample  matrix  at  the  output  of 
the  receive  array  has  the  form 

X  =  X«(0,)aRl(0,)a^(0I)S,+Xa(0e)aRl(0c)aL(0c)Sc+N  (1) 

t  C 

where  subscripts  t  and  c  index  into  the  set  of  targets  and  clutter  patches,  respectively.  Here,  a(0, ) 
represents  a  complex  scale  factor  associated  with  target  /,  while  i\Rx  (0, )  denotes  the  receive  array’s 
Nr  x  1  response  vector  for  target  t.  The  NtxR  matrix  S, ,  on  the  other  hand,  describes  how  the  various 
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waveforms  propagate  to,  and  are  reflected  by,  target  t.  That  is,  the  «lh  column  of  S,  contains  the 
convolution  of  sn  with  the  target’s  impulse  response.  Furthermore,  since  R»L,  S,  will  also  contain 
many  zeros  to  account  for  propagation  delays  to/from  the  target.  Variables  ) ,  aRi  (©,  ) ,  arr  (0r) 
and  St.  are  defined  similarly  with  respect  to  the  clh  clutter  patch.  Lastly,  the  NrxR  matrix  N  represents 
noise  as  well  as  other  interference. 


After  being  received,  the  sample  matrix  is  processed  to  extract  target  detections  and/or  to  form 
images.  Target  detection,  for  example,  can  be  performed  by  passing  the  received  samples  through  a  bank 
of  matched  filters.  Often,  this  filtering  is  implemented  in  stages,  wherein  the  initial  stage  applies  a  set  of 
“waveform”  matched  filters  to  form  V,  x  R'  matrices, 


Y_  = 


x  *lr 

Am  "l 


*hl, 


,  m  =  [,...,Nr 


(2) 


where  \rm  denotes  the  mlh  row  of  X,  *  denotes  convolution,  and  hn  denotes  the  Lxl  matched  filter  for 
signal  sn .  After  matched  filtering,  the  rows  of  the  NtNr  x R'  matrix  Y  =  T  Y,r  Y,y_  j  represent 
NtNr  channels  of  data.  It  is  thus  possible  to  coherently  combine  these  channels, 


z  =  w"Y  (3) 

to  achieve  some  desired  goal  (note:  incoherent  combining  of  channels  has  also  been  used  [9,  17]; 
however,  this  approach  will  not  be  discussed  further  here).  Often,  w  is  to  chosen  to  maximize  SINR,  i.e., 


Kvwv 


max- 


w 


"Rw 


where  v, ,  is  the  target  response  vector.  The  solution  to  this  problem  is  known  to  be  w  =  yR  v(f,( . 
R  =  £{  Y0  Y0" }  =  Rc  +  Rn  ,  Y0  are  target-free  snapshots  formed  in  a  manner  similar  to  Y,  Rc 
clutter  covariance  matrix,  RN  is  a  noise  covariance  matrix,  and  y  is  a  nonzero  complex  constant. 


(4) 

where 
is  the 


3.0  Doppler-Division  Multiple  Access  Waveforms 

Adaptive  MIMO  radar  performance  depends  on  the  availability  of  “orthogonal”  transmit  signals. 
To  date,  much  research  has  been  conducted  to  identify  waveform  sets  that  are  approximately  orthogonal. 
Toward  that  end,  one  approach  that  is  often  discussed  is  Frequency  Division  Multiple  Access,  or  FDMA. 

In  the  MIMO  radar  literature,  the  term  FDMA  is  used,  generically,  to  refer  to  sets  of  waveforms 
occupying  different  frequencies  at  the  same  time.  Mathematically,  the  //"’  FDMA  waveform  can  be 
written  as: 

sn(t)  =  s(t)eJ2*f’'  0  <t<T.  (5) 

for  n  =  \,...,Nt,  where  s(t)  is  a  common  radar  waveform  of  bandwidth  B,  and  fn  is  the  carrier 
frequency  for  waveform  n.  Usually,  the  carrier  frequencies  are  stepped  linearly,  e.g., 

fn  =fo  +  {n  "OAf  (g) 

where  f0  is  the  base  carrier,  and  Af  is  the  frequency  step  size.  In  the  context  of  MIMO  radar,  such 
waveforms  are  considered  quasi-orthogonal  as  long  as: 

i.e.,  Sj  j  ( t,rc )  must  be  negligibly  small  (with  respect  to  noise)  whenever  i  ^  j  . 
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On  transmit,  FDMA  M1MO  radars  emit  all  Nt  frequency  offset  waveforms  simultaneously. 
Then,  each  of  the  Nr  received  signals  is  processed  via  a  bank  of  matched  filters  -  one  filter  for  each 
waveform.  Each  filter,  in  essence,  demodulates  a  single  transmit  signal  and  compresses  it  in  range. 

There  are  many  types  of  FDMA  waveforms.  For  simplicity,  we  classify  these  into  three  groups 
(“noninterleaved,”  “coarsely  interleaved,”  and  “finely  interleaved”)  according  to  the  separation  between 
their  spectral  bands,  as  illustrated  in  Figure  I  (for  the  case  of  N,  =  2). 


(a) 


Frequency 

(b) 


Figure  1.  Examples  of  FDMA  waveforms. 

(a)  Non-interleaved  waveforms,  (b)  Coarsely  interleaved  waveforms, 
(c)  Finely  interleaved  waveforms. 


Noninterleaved  FDMA  waveforms  (see  Figure  la)  employ  frequency  offsets  that  are  larger  than 
the  signal  bandwidth,  i.e.,  |A/|  >  B ,  to  separate  the  transmitted  waveforms  during  each  transmit 

pulsewidth  (i.e.,  T  =  r).  Interleaved  FDMA  waveforms,  on  the  other  hand,  employ  smaller  frequency 
offsets  within  an  overlapping  frequency  band  (see  Figure  lb-c).  In  either  case,  frequency  offsets  must  be 
large  enough  to  isolate  the  Doppler  bandwidth  of  the  signals  of  interest  -  otherwise  ambiguities  will 
occur. 


By  design,  all  three  of  these  basic  FDMA  waveform  types  are  “orthogonal”  is  the  sense  of  (7). 
Unfortunately,  a  MISO  Cancellation  Ratio  (MCR)  analysis  [18]  reveals  a  significant  limitation  for 
FDMA.  In  particular,  the  adaptive  cancellation  performance,  as  measured  by  MCR,  degrades  quickly  as 
the  frequency  offset  grows,  according  to: 


MCR = 1 - 


sin(/rA//fi) 

xAf/B 


(8) 


In  fact,  it  is  clear  from  (8)  that  the  only  way  to  achieve  good  MCR  is  to  interleave  the  FDMA  waveform 
bands,  with  fine  interleaving  preferred  over  coarse  interleaving. 


One  approach  toward  such  fine  interleaving  is  to  use 

vsf./N,  (9) 

Then,  each  of  the  /V,  transmit  waveforms  will  occupy  a  “slice”  of  the  radar’s  unambiguous  Doppler 
band),  as  shown  in  Figure  1c.  On  receive,  the  MIMO  processor  can  then  use  conventional  Dopplers 
filters  to  isolation  each  transmit  channel.  Such  waveforms  are  called  Doppler  Division  Multiple  Access 
(DDMA). 
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In  practice,  the  phase  modulation  defined  by  (5)-(6)  and  (9)  is  approximately  constant  over  short 
time  intervals  such  as  the  radar  pulsewidth  T.  Hence,  the  DDMA  frequency  modulation  can  be 
approximated  via  a  piecewise  constant  phase  modulation.  That  is,  the  n'h  frequency  shifted  waveform  is 
created  by  applying  a  phase  modulation  to  the  common  waveform 

sn(t)  =  s{t)eK{p)  ~s{t)ei2nf-'  0 <t<T  (10) 

where  un  ( p)  =  2n  fn p / fr  denotes  the  phase  value  applied  to  transmitter  n  during  PRI  p  =  \_fr  f  J .  This 

approximation  enables  simplified  array  hardware,  since  each  transmit  channel  requires  only  a  phase 
shifters  (rather  than  a  means  of  independently  generating  frequency  offset  waveforms).  More 
importantly,  by  the  analysis  of  [18J,  excellent  clutter  cancellation  can  be  achieved. 

In  the  discussion  that  follows,  we  focus  on  DDMA  waveforms  as  defined  by  (6)  and  (9)  -  (10). 
However,  many  of  the  issues  we  describe  will  apply  to  systems  employing  true  (i.e.,  continuously 
modulated)  DDMA  as  well. 

3.1.  DDMA  Ambiguities 

To  receive  signals  (including  both  targets  and  clutter)  unambiguously,  the  DDMA  PRF  must 

satisfy 

fr/N,>Af>BD.  (11) 

For  large  Nt  and  Bu ,  such  DDMA  waveforms  will  have  relatively  high  PRFs.  In  Ground  Moving 
Target  Indication  (GMTI)  radar,  for  example,  these  PRFs  (/r  >  N,BD)  will  be  up  to  Nt  times  higher  than 
“normal”  SIMO  GMTI  radar  PRFs  ( fr  >  BD ).  As  a  result,  there  can  be  a  larger  number  of  both  “range 
ambiguities”  and  range  eclipsed  regions' . 

To  address  the  problem,  we  could  lower  the  PRF  to: 

fr/Nt>Af>Bc  (12) 

(where  Bc  is  the  Doppler  bandwidth  associated  with  clutter),  resulting  in  PRFs  that  (1)  are  high  enough 
so  as  to  unambiguously  isolate  the  clutter  from  each  transmit/receive  signal  path,  but  (2)  allow  high 
speed  targets  to  become  Doppler  ambiguous,  hi  this  case,  a  high  speed  target  signal  received  via  one 
transmit/receive  signal  path  will  look  the  same  as  a  low  speed  target  signal  (or  clutter  signal)  received  via 
a  different  transmit/receive  signal  path.  When  this  happens,  high  speed  targets  will  “fold  over”  onto  low 
speed  targets.  Moreover,  some  high  speed  targets  will  alias  into  the  blind  zones  associated  with  the  low- 
speed  clutter. 

Lowering  the  PRF  further,  to  fr/Nt  >Af>eBc  with  e<\,  leads  to  additional  problems.  Such 
low  PRFs  increase  the  clutter  rank  and  require  additional  degrees  of  freedom  for  clutter  suppression. 
Hence,  most  systems  will  choose  to  operate  in  the  PRF  regime  dictated  by  (12)  -  which  results  in 
Doppler  ambiguous  high  speed  targets. 

3.2.  Impact  of  Ambiguities 

The  performance  of  DDMA  using  (12)  is  illustrated  in  a.  Here,  a  sidelooking  airborne  GMTI 
radar  is  simulated.  The  radar  uses  Nt  =16  transmit  elements,  arranged  linearly  with  an  inter-element 
spacing  of  dF  =0.19  m,  and  N r  =  1  receiver  element.  The  element  gain  is  13  dBi  and  the  wavelength  is 
X  =  0. 1 265  m.  The  radar  platform  is  moving  at  v  =  27.5  m/s  at  an  altitude  of  10000  ft. 


f  Range  ambiguities  and  eclipsing  could  be  reduced  by  re-designing  the  radar  to  smaller  /V,  or  II n  .  For  the  remainder  of  Uiis  paper,  however, 
N ,  and  11 D  are  assumed  fixed. 
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According  to  the  theory  of  airborne  adaptive  radar  [19],  this  radar’s  mainlobe  clutter  should  have 
a  Doppler  spread  of  ±  2vcos(£o)sin(0fl)//l .  In  our  simulation,  the  ground  was  illuminated  at  a  slant 
range  of  6  km  (corresponding  to  a  depression  angle  of  £D  =  30° )  and,  due  to  element  pattern  effects,  the 
radar’s  mainlobe  clutter  power  fell  below  noise  at  an  azimuth  offset  of  0B=34°.  Consequently,  the 
Doppler  spread  is  ±2vcos(£D)sin(i9B)//l  =  ±  210.5  Hz,  corresponding  to  a  Doppler  bandwidth  of 
Bc  =  421  Hz.  Had  our  system  used  a  single  transmit  phase  center  (i.e.,  SIMO  radar  with  Nt=\),  all  this 
clutter  would  appear  in  a  single  Doppler  band  (spread  over  +210.5  Hz).  With  N ,  =16  transmitters, 
however,  the  clutter  spectrum  contains  multiple  clutter  bands,  as  shown  in  a.  In  the  MIMO  case,  each  of 
these  Nt-\  “extra”  clutter  bands  is  associated  with  one  of  our  Nt-l  Doppler-offset  transmit 
waveforms.  In  particular,  the  /i,h  transmitter’s  clutter  is  shifted  by  (n  -  1)  A/  Hz.  In  our  simulation,  we 
used  a  PRF  of  /r=  6745  Hz  and  a  frequency  offset  of  Af  =  fj Nt  =421.6  Hz,  yielding  a  shift  of 
(n -l)A/  =  421. 6n -421.6  Hz.  This  relationship  (i.e.,  the  mapping  of  transmitters  to  clutter  bands)  is 
made  explicit  via  the  annotations  at  the  top  of  a. 

For  GMTI  radar,  the  multi-band  clutter  spectrum  (depicted  in  a)  must  be  cancelled  to  detect 
moving  targets.  Fortunately,  since  our  frequency  offset  (A/  =421. 6 Hz)  is  larger  than  the  total  clutter 
bandwidth  (Bc  =421  Hz),  the  clutter  return  from  each  transmitter  is  well-isolated  in  the  Doppler 
domain.  This  isolation  allows  the  MIMO  signal  processor  to  unambiguously  isolate  the  clutter  due  to 
each  transmit  phase  center.  The  MIMO  receiver  can  then  use  these  degrees  of  freedom  to  cancel  the 
clutter. 


Fast  targets,  however,  will  remain  Doppler  ambiguous  and  will  thus  alias  back  into  the  endo- 
clutter  region,  b  -  c  depicts  the  resulting  SINR  losses  after  adaptive  clutter  cancellation- .  Note  that  the 
maximum  possible  clutter  Doppler  shift,  2v  cos  (£",,)  sin  (<9B)//l  =  210.5  Hz,  occurs  when  clutter  is 
received  from  the  leading  edge  of  the  element  pattern.  Likewise,  the  minimum  Doppler  shift, 
-2v  cos  (  (fD  )  sin  (#B)//t  =  -  210.5  Hz,  occurs  when  clutter  is  received  from  the  aft  edge  of  the  element 
pattern.  At  all  other  angles  in  between,  the  clutter  Doppler  varies  according  to 

2 v cos  ( £. ) sin  ( dc  )/X ,  (13) 

forming  a  “ridge”  when  viewed  in  the  Doppler-spatial  frequency  plane.  Now,  had  our  system  used  just  a 
single  transmit  phase  center  (N,  =1),  (13)  dictates  the  existence  of  a  single  clutter  ridge.  SINR  losses 
would  be  high  along  this  ridge,  but  close  to  zero  elsewhere.  With  Nt  =16  Doppler-offset  transmit 
waveforms,  however,  the  aggregate  frequency  shift  associated  with  transmitter  n  and  clutter  patch  c  will 
be 


(n  - 1)  Af  +  2vcos  ( £ )  sin  (0c  )/A  ( ,  4) 

Consequently,  there  will  be  N ,  clutter  ridges.  This  results  in  N,- 1  additional  ridges  of  high  SINR  loss, 
as  shown  in  b§.  This  also  means  there  will  be  IV,—  1  new  “blind  speeds”  caused  by  the  DDMA 
waveform  at  any  given  target  angle,  c  shows  a  cut  through  the  SINR  loss  plot  at  array  broadside.  There, 
we  see  Nt  narrow  Doppler  regions  (i.e.  target  speeds)  having  high  loss  due  to  the  clutter  nulls  formed 
during  STAP  processing.  Comparable  low-PRF  SIMO  radars,  in  contrast,  would  have  only  a  single  blind 
speed. 


*  To  make  these  plots,  the  N,Nr  MIMO  channels  were  adaptively  weighted  and  combined  using  Space-Time  Adaptive  Processing  (STAP). 
STAP  is  a  variation  on  the  SINR  maximization  process  of  (3)  -  (4)  used  to  adaptively  cancel  clutter  [WARD].  In  a  STAP  system,  adaptive 
clutter  nulls  are  automatically  adjusted,  as  a  function  of  Doppler,  to  track  clutter.  Losses  will  occur  in  regions  that  are  close  to  the  clutter 
manifold.  These  “SINR  Losses,"  which  are  depicted  in  b  -  c,  are  computed  as  the  ratio  of  the  STAP  processed  output  to  the  ideal,  nonadaptivcly 
processed  output  in  an  interference  free  environment.  Hence,  a  loss  of  0  dB  is  ideal,  while  large  negative  values  indicate  regions  of  high  loss. 

5  Note  that  b  plots  loss  as  a  function  of  Doppler  frequency  and  normalized  spatial  frequency,  with  normalized  spatial  frequency  defined  as 
'/ccos(£)sin(<9,)/'*  >■ 
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Figure  2.  GMTI  radar  using  DDMA  waveform  ( Nt  =  16,  N r  =  1). 

(a)  Clutter  spectrum,  (b)  SINR  Loss  after  STAP  processing  (we  used  the  PRI 
staggered  post-Doppler  STAP  algorithm  of  [19]),  (c)  SINR  loss  at  broadside 


3.3.  Analysis  of  Ambiguities 

The  Nt- 1  extra  clutter  ridges,  discussed  in  Section  3.2,  occur  because  fast  targets  (i.e.,  any 
target  t  having  Doppler  shift  d,e.  DtAST  ={d,\ |r/,|>  fr/2N,])  alias  in  Doppler.  In  some  cases,  the 
aliased  Doppler  will  be  on  (or  near)  the  clutter  ridge.  When  this  happens,  the  STAP  weights  used  to 
receive  a  target  t  (as  defined  by  parameters  0,  =  { 0, ,  £ ,  dt  e  DFASJ ,  z( } )  are  nearly  the  same  as  the 
weights  used  to  receive  a  slow  target  t'  (as  defined  by  0'  =ldl,£l,d' ,z,  where 
r/,  e  oslow  s{d':\d'\<fr/2N,} ) -  resulting  in  a  loss  when  the  later  is  close  to  the  true  clutter  ridge  (aka 
blind  speed). 

In  SIMO  radar,  this  sort  of  aliasing  occurs  when  the  target  Doppler  is  \d\>  fr/2  Such  targets 
are  called  “Doppler  ambiguous”.  In  DDMA  MIMO  radar,  Doppler  aliasing  occurs  at  slower  target 
speeds,  i.e.,  whenever  |r/,|>  fr/2Nt.  This  “early-onset”  Doppler  aliasing  is  a  result  of  the  DDMA 
MIMO  processing. 

There  is  one  important  difference,  however,  between  DDMA  Doppler  aliasing  and  the  SIMO 
Doppler  aliasing.  Specifically,  in  DDMA,  fast  target  echoes  from  the  /ilh  transmit  signal  will  be  received 
and  processed  in  such  a  way  as  to  become  Doppler  ambiguous  with  slow  targets 
(d'e  Dsww  ={d':  \d'\  <  fr  j2Nt } )  from  a  different  MIMO  transmitter.  This  relationship  is  made  explicit 
by  writing  the  transmit  array  response  vector  as  a  function  of  0'  ={0l,^l,c l',z,} ,  which  includes  several 
additional  target  parameters  such  as  the  target  Doppler  d'  and  target  range  zone  z,  ■  For  DDMA, 
aTx[0t’Ci’d,e  DtAsr  ,z,)  is  a  permutation  of  an (#,,£, rf/e  DSLOW,z,).  This  property  results  in  a  unique 
kind  of  ambiguity,  as  illustrated  in  Figure  3. 

The  exact  way  in  which  elements  of  the  ar<  ((?,,£, r/'e  Dslow ,z,)  are  permuted  to  form  the 
elements  of  a Tx[9r£t,dte  DfAST ,z,)  depends  on  the  order  that  the  DDMA  frequencies  are  mapped 
to  transmitters.  The  standard  linear  mapping,  (6),  results  in  a  cyclic  permutation.  That  is, 
ay> {9t’C,’d,e  DIASr ,z,)  is  formed  by  circular  shifting  (or  rotating)  of  art  ( #, , £ , d' e  DSI0W , z( ) .  The 
size  of  this  shift,  Nshi/I ,  depends  on  the  size  of  dt  e  DfAST .  If  |r/(|  is  large,  the  shift  will  be  large;  if  |r/,|  is 
small,  the  shift  will  be  small.  Quantitatively,  we  may  define: 

Nm/,  =  s‘gn  (<'.)L(KI+/./2".)/(/./'v.)J  os) 

and 
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(16) 


where  [  J  denotes  the  floor  operation,  sign(x)  =  1  when  .v  >  0 ,  and  sign(x)  =  -l  if  .v < 0 .  Then, 
a  Tx(e„^d,eDFAS\Zl)  is  formed  by  circular  shifting  Nshjft  times,  i.e., 


Qhtfi  Ku  M 


(17) 


For  uniform  linear  arrays,  this  permutation  causes  strong  “ambiguities.”  For  uniform  linear  arrays,  the 
target  response  vector  when  d'e  d5LOW  takes  the  form: 


(18) 


where  t//  =  2dHcos(^l)s\n(0l)/ A ,  Note  that  the  elements  of  (18)  have  linear  phase  slope.  For  fast 
targets,  (18)  is  circularly  shifted  Nshifl  times.  Nonetheless,  the  shifted  vector  still  has  the  same  linear 
phase  slope,  except  for  a  phase  discontinuity  at  element  NMfl .  Depending  on  the  location  of  this  phase 
discontinuity,  the  resulting  array  response  can  vary  from  a  “sum”  to  “difference”  pattern  when  computed 
over  the  non-rotated  manifold.  Regardless  of  discontinuity  location,  the  resulting  response  vector  is 
highly  correlated  with  slow-speed  targets  near  the  same  angle,  resulting  in  a  “near”  ambiguity.  At  certain 
target  Dopplers,  moreover,  this  linear  phase  progression  will  be  nearly  the  same  as  that  of  clutter.  As 
such,  clutter  nulls  can  appear  ambiguously  in  Doppler. 

3.4.  Avoiding  Ambiguities  with  Normal  DDMA 

Nonuniformly  sampled  linear  (or  planar)  arrays,  of  course,  do  not  have  response  vectors  of  the 
form  (18).  Nonplanar  arrays  are  also  unlikely  to  have  response  vectors  of  the  form  of  (18). 
Consequently,  such  array  topologies  provide  one  approach  toward  mitigating  the  Doppler  ambiguity 
problem  described  in  Section  3.3. 

Unfortunately,  most  radars  use  uniform  linear  (or  uniform  planar)  array  topologies.  Such 
topologies  dominate  the  marketplace  due  low  manufacturing  costs  and  many  other  practical  reasons. 
Hence,  other  approaches  toward  mitigating  the  Doppler  ambiguity  problem  are  needed.  Such  techniques 
are  the  focus  of  the  next  section. 
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Figure  3.  Transmitter-dependent  phase  values  in  art(0, )  for  linearly  mapped 

DDMA  waveform  (6). 


4.0  Dithered  DDMA 

To  eliminate  the  DDMA  blind  speeds  described  in  Section  3.3,  we  developed  two  modified 
DDMA  techniques.  Collectively,  we  call  these  “dithered  DDMA”  techniques.  The  first  technique  works 
by  altering  the  mapping  between  DDMA  transmit  frequencies  and  transmit  antennas.  The  second 
technique,  in  contrast,  adds  a  phase  offset  to  each  of  the  transmitted  waveforms.  In  both  cases,  these 
techniques  work  by  altering  the  array  response  vectors  for  fast  targets,  thereby  suppressing  ambiguities. 
Better  still,  these  techniques  are  complementary,  meaning  they  can  be  used  together  for  even  greater 
benefit. 


4.1.  Frequency-Dithered  DDMA 

In  a  frequency-dithered  DDMA  system,  we  alter  the  mapping  of  DDMA  frequencies  to  transmit 
antennas.  This,  in  turn,  modifies  the  permutation  so  that  it  can  no  longer  be  represented  as  a  series  of 
circular  shifts,  i.e.: 

.  .  Permutation  *  Shift  by 

»t -*  *0-  09) 
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The  goal  is  to  choose  a  mapping  that  alters  a,.t(0, )  in  some  favorable  way.  The  mapping  could,  for 

example,  be  pseudo-random  (which  provides  good  performance  over  all  angles).  The  mapping  could  also 
be  optimized  to  favor  specific  target  angles  and  speeds.” 

The  impact  of  pseudo-random  frequency-dithering  is  clearly  illustrated  in  Figure  4.  Here,  we 
plot  the  SINR  loss  for  an  airborne  GMTI  radar  having  the  same  parameters  as  in  .  Away  from  the 
broadside  region  (i.e.,  toward  the  left  or  right  sides  of  Figure  4a),  the  impact  of  frequency-dithering  is 
particularly  evident.  In  these  regions,  the  array  response  vector’s  phase  slope, 
^  =  2r/fc.cos(£)sin(0()/,2. ,  will  be  large  because  |^|»0.  Random  permutations  of 

a7r DSI0W , z, )  thus  result  in  response  vectors,  a ^ (#,,£, r/f  e  DFAST ,z,) ,  that  generally  have  a 
small  projection  onto  their  unpermuted  counterparts,  aTx  (#f,£,,r/,'e  Ds,ow  ,z,) .  Consequently,  this 
allows  good  nulling  of  clutter  at  the  same  angle  as  the  fast  target. 

Near  the  broadside  region,  however,  an  ,  £  ,rf'e  dslow  ,  z, )  has  a  small  phase  slope  \if/\. 
Consequently,  random  permutations  of  will  result  in  response  vectors, 

aTx{fy’Ci’die  DFAST>z,),  that  are  much  closer  to  their  unpermuted  counterparts.  In  the  extreme  (i.e.,  at 
array  broadside,  where  ari(^,C,,rf'G  Ds'w,z()  =  [l  •••  l]r  and  ^  =  0),  this  distance  becomes  zero; 
hence,  permuting  a7 rx(Ot,£t,d’te.  DSLOW ,zt)  has  no  effect  on  it.  Therefore,  it  can  be  difficult  (or 

impossible)  to  detect  fast  targets  near-broadside  region.  This  point  is  illustrated  in  Figure  4b,  which 
shows  the  SINR  loss  at  broadside.  Note  that  Nt  blind  speeds  are  present,  despite  frequency-dithering. 
Another  technique  must  be  used  to  suppress  blind  speeds  in  the  broadside  region. 


Figure  4.  GMTI  radar  using  Frequency-Ditliered  DDMA  waveform, 
(a)  SINR  Loss,  (b)  SINR  loss  at  broadside 


”  For  example,  the  mapping  { /, ,  fs, ,  /3 .  /v_. { '/a  i . 7 \2 , v:v3, 7\-4, . . . }  shifts  targets  in  Doppler  zone  2  so  that  they  compete 
with  clutter  from  angle  -e, .  This  sort  of  mapping  might  be  used  when  the  transmit  antenna  has  low  gain  toward  -o, 
(because  targets  in  Doppler  zone  2  will  compete  with  clutter  that  is  weak  and  unlikely  to  cause  a  blind  speed).  There  are 
many  other  similar  specialized  mappings  that  might  be  considered. 
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4.2.  Phase-Dithered  DDMA 


Near  the  broadside  region,  frequency-dithered  DDMA  fails  to  mitigate  the  DDMA  blind  speeds 
discussed  in  Section  3.3.  Fortunately,  a  second  technique,  called  “phase-dithered  DDMA,”  works  well  in 
this  region. 

In  a  phase-dithered  DDMA  system,  we  add  a  phase  offset  to  each  of  the  DDMA  transmit  signals. 
That  is,  we  transmit 


»„(<)  =  *(  (20) 

Note  that  the  ii  transmitter  has  been  assigned  a  phase  value  vn . 

On  receive,  the  MIMO  processor  applies  matched  filters  for  both  the  frequency  and  the  phase 
of  the  transmitted  signal.  For  d ,  e  dslow  >  this  process  removes  the  transmit  phase  vn  (i.e.,  it  multiplies 
the  received  signal  by  e~jv"  ).  For  dt  e  DFAST ,  however,  the  target  signal  will  appear  at  the  output  of  the 
“wrong”  matched  filter  (e.g.,  it  appears  in  matched  filter  n  +  Nshjfl ,  assuming  a  linear  DDMA  mapping  is 

used  as  in  (6))  due  to  the  Doppler  ambiguity  effect.  As  discussed  in  Sections  4.0  and  4.1  ,  this  will  cause 
the  target  to  have  a  permuted  transmit  array  response  vector.  Likewise,  it  will  have  the  “wrong”  set  of 

phases  removed  (i.e.,  it  will  be  multiplied  by  e  JV"‘Sik#  (  resulting  in  a  phase  residue  of  Vn-Vn+Nur ). 

Choosing  vn  to  be  uniformly  random  on  [0, 2/r]  makes  these  phase  residuals  random,  providing  robust 
performance  over  all  angles.  Other  choices  of  vn  can  be  used  to  optimize  performance  over  smaller 
angular  sectors. 


Interestingly,  this  idea  of  inserting  a  spatially  varying  phase  dither  and  then  later  removing  it, 
with  the  net  effect  of  shifting  undesired  signals  away  from  the  desired  signal  subspace,  is  not  new.  This 
same  idea  has  been  used,  previously,  to  mitigate  nonlinear  distortion  effects  originating  in  the  radar 
electronics  [20].  For  this  work,  we  merely  adapted  the  concept  to  MEMO  radar. 


Figure  5  illustrates  the  impact  of  phase-dithering  applied  to  DDMA.  Here,  the  N,  phase  values, 
{u,,} ,  were  taken  to  be  uniformly  random  on  [0,2/r] .  Pseudo-random  frequency-dithering  was  also  used. 

Note  that  the  losses  due  to  ambiguous  clutter  near  boresight  (Figure  5b)  have  been  replaced  with  a  much 
smaller  loss  spread  over  the  entire  Doppler  space. 


Normalized  Spatial  Frequency 


(a)  (b) 

Figure  5.  GMTI  radar  using  Phase-Dithered  DDMA  waveform, 
(a)  SINR  Loss,  (b)  SINK  loss  at  broadside 
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4.3.  Range  Ambiguities 

In  pulse-Doppler  radar,  range  ambiguities  occur  whenever  there  are  detectable  echoes  from 
targets  (or  clutter)  at  ranges  greater  than  c/2/r  .  Ranges  under  c/2/r  are  said  to  be  unambiguous, 
whereas  ranges  z,c/ 2 fr<pt<(z,  + l)c/2 fr  are  said  to  lie  in  the  z,  -th  ambiguous  range  zone. 

In  DDMA  MIMO  radar,  range  ambiguous  response  vectors  art(0;, £,</,, z, )  will  differ  from 
unambiguous  responses  (0(,  £,</,,()) ,  even  if  both  are  at  the  same  angle  and  Doppler.  To  detect  range 
ambiguous  targets,  then,  the  correct  steering  vector  should  be  used. 

For  linearly  mapped  DDMA  (as  in  (6)),  the  //"'  element  of  aTx  (<?,,£, d, , z, )  is  shifted  by 
as  ii|ustrated  in  Figure  3.  This  shift  has  an  interesting  effect  on  the  response  vector 
aTx(&t’Ci’di>z,)  ■  In  particular,  it  modulates  art(<?(,£,r/(,z() ,  causing  it’s  pattern  to  points  toward  a 
shifted  angle.  This  results  in  an  ambiguity  between  targets  (or  clutter)  in  zone  z,  and  targets  (or  clutter) 
in  zone  0  at  some  shifted  angle.  The  size  of  this  angular  shift  increases  with  z, ,  meaning  that  range 
zones  that  are  “close”  will  have  small  angular  shifts,  while  range  zones  that  are  far  apart  have  a  large 
angular  shifts.  This  could  be  an  issue;  we  are  often  interested  in  signals  received  from  a  small  number  of 
range  zones  that  are  relatively  close  together.  In  such  cases,  the  small  angular  shifts  of  the  response 
pattern  could  result  in  clutter  (or  targets)  from  nearby  angles  competing  with  the  signal  of  interest 

Frequency  dithered  DDMA,  however,  solves  this  problem.  By  altering  the  mapping  of 
frequencies  to  transmitters,  we  modify  the  modulation  applied  to  range-ambiguous  response  vectors. 
With  certain  frequency  dithers,  for  example,  we  can  increase  the  size  of  the  angular  shift  applied  to 
nearby  range  zones77.  This  type  of  dithering  might  be  used  to  better  separate  signals  and  clutter  in  nearby 
range  ambiguities. 

Likewise,  we  could  pseudo-randomly  map  frequencies  to  transmitters.  This  type  of  dithering 
results  in  modulation  sequences  that  do  not  shift  the  angle  response.  Instead,  they  impart  unique  phase 
signatures  onto  each  range  zone’s  response  vectors.  The  response  vectors  from  different  range  regions 
will  thus  have  low  cross-correlation.  Furthermore,  these  response  vectors  can  be  used  as  steering  vectors 
by  and  adaptive  clutter  cancellation  algorithm  (like  STAP)  -  further  suppressing  range  ambiguities 
(much  like  the  way  Doppler  ambiguities  were  suppressed  in  Section  4.1). 

5.0  Summary  and  Conclusions 

In  this  paper,  we  examined  one  promising  type  of  FDMA  waveform,  in  which  disjoint  Doppler 
frequency  bands  were  used  to  isolate  the  transmitted  signals.  In  particular,  these  “Doppler  Division 
Multiple  Access”  waveforms  were  shown  to  support  high  MCR  levels.  However,  these  waveforms  also 
suffered  from  clutter  ambiguities  and  blind  speeds.  The  source  of  these  ambiguities  was  investigated. 
Their  origin  was  shown  to  be  the  use  of  linear  frequency  shifts  across  a  uniform  linear  antenna  aperture. 

To  address  this  problem,  we  proposed  a  modified  DDMA  waveform  approach.  In  part,  this 
approach  leveraged  prior  work  relating  to  mitigation  of  array  distortion  effects.  In  simulations,  this 
modified  DDMA  approach  was  able  to  eliminate  blind  speeds  and/or  shift  them  to  regions  of  the  Doppler 
space  where  they  are  more  tolerable. 
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